Eccentricity fluctuation effects on elliptic flow in relativistic heavy ion collisions 
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We study effects of eccentricity fluctuations on the elliptic flow coefficient V2 at mid-rapidity in 
both Au+Au and Cu+Cu collisions at ^/sjvjv = 200 GeV by using a hybrid model that combines 
ideal hydrodynamics for space-time evolution of the quark gluon plasma phase and a hadronic 
transport model for the hadronic matter. For initial conditions in hydrodynamic simulations, both 
the Glauber model and the color glass condensate model are employed to demonstrate the effect 
of initial eccentricity fluctuations originating from the nucleon position inside a colliding nucleus. 
The effect of eccentricity fluctuations is modest in semicentral Au-|-Au collisions, but significantly 
enhances V2 in Cu+Cu collisions. 
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I. INTRODUCTION 

One of the major discoveries at Relativistic Heavy 
Ion Collider (RHIC) in Brookhaven National Laboratory 
(BNL) is that azimuthally anisotropic flow (the so-called 
elliptic flow) [l| is found to be as large as an ideal hy- 
drodynamic prediction for the first time in relativistic 
heavy ion collisions [1, d, [I]- Whether local thcrmal- 
ization, which is demanded in application of hydrody- 
namics, is reached is not known a priori in relativistic 
heavy ion collisions. Therefore the discovery indicates 
that the heavy ion reactions at relativistic energies pro- 
vide a unique opportunity to investigate the high tem- 
perature QCD matter in equilibrium. The discovery also 
indicates a new possibility to obtain the transport prop- 
erties of the QCD matter under extreme conditions. Ef- 
fects of viscosity have already been taken into account in 
several hydrodynamic simulations d, @, 0, H, Q and also 
specifically in distortion of distribution functions [HI 
to calculate elliptic flow coefficients. 

Systematic studies showed, however, that the reason- 
able agreement between results from ideal hydrodynam- 
ics and elliptic flow data has been achieved only by a 
particular combination of dynamical modeling, namely 
initial conditions from the Glauber model, the perfect 
fluid quark gluon plasma (QGP) core, and the dissipa- 
tive hadronic corona For example, when one takes 
initial conditions from the color glass condensate (CGC) 
model instead of the conventional Glauber model ones, 
the result overshoots the elliptic flow data due to ini- 
tial eccentricity larger than that from the Glauber model 
p^ . Even within the Glauber model initializations, the 
agreement between hydrodynamic results and the data, 
in particular centrality dependence of the elliptic flow, 
is not perfect, which would be due to an absence of ini- 
tial fluctuation effects 0, [H, [13, [H, [H, M, [Ml, [H, [H] . 
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Therefore, further investigation is indispensable toward 
better understanding of the elliptic flow data and, in 
turn, understanding of transport properties of the QGP. 

In this paper, we follow up with the previous study 
[T3I . [2^ 1 based on an ideal hydrodynamic description 
of the QGP fluid and a microscopic description of the 
hadronic gas by taking into account the initial eccen- 
tricity fluctuation effects. It was found that large ini- 
tial eccentricity in some particular modeling of the CGC 
picture is attributed to improper treatment of nuclear 
edge regions to define the saturation scale [20| . Improved 
treatment of the edge leads to slight reduction of eccen- 
tricity [lOl . So it would be interesting to see whether the 
improved CGC model ends up with reproduction of ellip- 
tic flow data as well as how initial fluctuation of eccen- 
tricity affects elliptic flow coefficients in both the Glauber 
model and the CGC model. 

This paper is organized as follows. In Sec. HH we first 
overview the dynamical modeling of relativistic heavy ion 
collisions based on a hybrid (hydro-|-cascade) approach. 
We discuss how to implement the eccentricity fiuctuation 
in the hydrodynamic initial conditions. For the initial 
transverse profiles, we employ the Glauber model and 
the CGC model and compare the eccentricities with each 
other. In Sec. IIIIl we investigate the effect of eccentric- 
ity fluctuation on the elliptic flow coefficients by using 
the hybrid dynamical model. Section |IV] is devoted to 
conclusion. 



II. THE MODEL 

A. Dynamical modeling of heavy ion collisions 

Our study is mainly based on a hybrid model which 
combines an ideal fluid dynamical description of the QGP 
stage with a realistic kinetic simulation of the hadronic 
stage [H [H [H [H, [13, [H, [ll]. Relativistic hydro- 
dynamics is the most relevant framework to understand 
the bulk and transport properties of the QGP since it 
directly connects the collective flow developed during 
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the QGP stage with its equation of state (EOS). It is 
based on the key assumption of local thermalization. 
Since this assumption breaks down during both the very 
anisotropic initial matter formation stage and the di- 
lute late hadronic rescattering stage, the hydrodynamic 
framework can be applied at best only during the inter- 
mediate period. To describe the breakdown of the hydro- 
dynamic description during the late hadronic stage due to 
expansion and dilution of the matter, one may have two 
options: One can cither impose a sudden transition from 
thermalized matter to non-interacting and free-streaming 
hadrons through the Coopcr-Fryc prescription [sot at a 
decoupling temperature Tdcc, or make a transition from a 
macroscopic hydrodynamic description to a microscopic 
kinetic description at a switching temperature Tgw We 
here use the first approach to fix the initial parameters by 
comparison with the multiplicity data. After all initial 
parameters are fixed, we carry out hybrid simulations to 
investigate the effect of initial eccentricity fluctuations as 
well as hadronic dissipation. 

For the space-time evolution of the perfect QGP fluid 
we solve numerically the equations of motion of ideal hy- 
drodynamics for a given initial state in three spatial di- 
mensions and in time [s^ : 

^^T^^•' = 0, (1) 
T^"' = {e + p)u''u'' - pg^"' . (2) 

Here e, p, and ti^ are energy density, pressure, and 
four-velocity of the fluid, respectively. We neglect the 
net baryon density due to its smallness at collider ener- 
gies Isil. We have solved Eq. ([1]) in (r, x, y, rjs) coordi- 
nates [32| where x± = {x,y) is a transverse coordinate, 
T = vt^-z^ and rjs — ^\n[{t+z)/ {t—z)] are longitudinal 
proper time and space-time rapidity, respectively, ade- 
quate for the description of collisions at ultra-relativistic 
energies. As will be discussed later in the next subsec- 
tion, we calculate the initial condition only at midra- 
pidity and neglect possible fluctuation and correlation 
effects in the forward/backward rapidity regions in this 
work. So we assume longitudinal boost invariance up to 
the longitudinal boundary of the mesh at 77^, = 5 in three 
dimensional grids and calculate the observables only at 
midrapidity. The solution for the positive rapidity region 
is properly reflected to the negative rapidity region. We 
have checked that there is neither rapidity dependence 
nor numerical artifacts from the finite volume in the lon- 
gitudinal direction on the observables at midrapidity. 

For the high temperature QGP phase [T > = 170 
MeV) we use the EOS of massless parton gas (w, d, s 
quarks and gluons) with a bag pressure B: 

p=i(e-4B) (3) 

The bag constant is tuned to = 247.19McV to ensure 
a first order phase transition to a hadron resonance gas 
at critical temperature Tc = 170 MeV. The hadron reso- 
nance gas model at T < Tc includes all hadrons up to the 



mass of the A(1232) resonance. Our hadron resonance 
gas EOS implements chemical freeze-out at Tchom = 7c = 
170 MeV [33|. This is achieved by introducing appropri- 
ate temperature-dependent chemical potentials fJ.i{T) for 
each hadronic species i [H, [H, [H, [H, [13, HI] . In our 
hybrid model simulations we switch from ideal hydro- 
dynamics to a hadronic cascade model at the switching 
temperature Tgv, = 169 MeV. The subsequent hadronic 
rescattering cascade is modeled by JAM [39|, initialized 
with hadrons distributed according to the hydrodynamic 
model output and calculated with the Cooper-Frye for- 
mula [lOl along the Tgw = 169 MeV hypersurface rejecting 
inward-going particles. 



B. Initial conditions with eccentricity fluctuation 
effects 

The space-time evolution of thermodynamic variables 
is described in the hydrodynamic framework. The con- 
cept of ensemble average in a sense of statistical mechan- 
ics is implicitly there to interpret the thermodynamic 
variables. So we identify the ensemble with a large num- 
ber of collision events and suppose the hydrodynamic so- 
lution represents an average behavior of the space-time 
evolution of the matter for a given centrality cut. This 
is contrary to an approach in Refs. (l9l . [40| in which hy- 
drodynamic equations with lumpy initial conditions are 
solved in an event-by-event basis. The smooth initial con- 
dition is, however, required in some practical reasons in 
our approach. Such a lumpy initial condition could gen- 
erate sizable numerical viscosity in hydrodynamic simu- 
lations if the mesh size is not sufficiently small. In the 
present study, the mesh size in the transverse direction 
is Ax = 0.3 (0.2) fm in Au-hAu (Cu-hCu) coUisions. The 
mesh size is of the same order of the transverse size of 
nucleons, so it would be hard to capture possible lumpy 
structures in the initial conditions. Moreover, one needs 
to perform a large number of simulations to gain suffi- 
cient statistics in final observables. Therefore we pursue 
a more conservative approach by initializing the distribu- 
tion of thermodynamic variables with smooth functions 
averaged over many events. 

For initial conditions in this study, we employ the 
Monte-Carlo version of both the Glauber model and the 
factorized Kharzcev-Lcvin-Nardi (fKLN) model [1^ to 
generate the initial distribution of entropy density in 
an event-by-event basis. One has extensively used the 
Monte Carlo Glauber model (MC-Glauber) to determine 
the centrality cut, the average numbers of participants 
and binary collisions for a given centrality, and so on 
[4l|. On the other hand, the fKLN model enables us to 
improve the treatment of entropy production processes. 
This model also gives a natural description near the edge 
regions compared to the ordinary KLN approach [i^ em- 
ployed by us previously [i^ . The Monte-Carlo version of 
fKLN, which we call MC-KLN, implements the fluctu- 
ations of gluon distribution due to the position of hard 



3 



sources (nucleons) in the transverse plane [20| • 

We first calculate a transverse entropy density profile 



sa{x±_) = s(t = To,x,y,ris = 0) 



(4) 



in each sample at an impact parameter b for a given cen- 
trality, where tq ~ 0.6 fm/c is the initial time for the 
hydrodynamical simulations. Then the variances of the 
profile are obtained from 



''xy = {xy) - (x){y). 



(5) 
(6) 
(7) 



Here (■ ■ • ) describes the average over transverse plane by 
weighting the entropy density in a single sample: 



/ £x± ■ ■ ■ so{x±) 
J d'^x±so{xj_) 



(8) 



The eccentricity with respect to reaction plane, the par- 
ticipant eccentricity, and the corresponding transverse 
areas can be defined [l^, respectively, as 



£rp — 

^part — 
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part T^sJ'^x'^y- CTxy 



(9) 

(10) 

(11) 
(12) 



It should be noted here that, in calculating eccentric- 
ity and transverse area, we use the entropy density pro- 
file instead of the distribution of participants. Thus the 
eccentricity depends on modeling of initialization, i.e., 
transverse profiles in hydrodynamic simulations. 

The impact parameter vector and the true reaction 
plane are not known experimentally. So one can set 
an apparent frame of created matter shifted by (x, y) = 
((x), (y)) and tilted by ^' from the true frame in the trans- 
verse plane (isj : 



tan 2* 



2a 



(13) 



xy 



The anisotropy of particle distribution could be corre- 
lated with the misidentified frame. To take account of 
this, we first shift the center-of-mass of the system to the 
origin in the calculation frame and then rotate the profile 
in the azimuthal direction by ^' to match the apparent re- 
action plane to the true reaction plane. We generate the 
next sample of an entropy density profile again as above 
and average the profiles over many samples. We repeat 
the above procedure for many samples until the initial 
distribution is smooth enough. The initial conditions ob- 
tained in this way contain the effects of eccentricity fluc- 
tuation even though the profile is smooth. In particular. 



even in case of vanishing impact parameter, eccentric- 
ity is finite due to its evcnt-by-event fluctuation. It is 
the particle distribution calculated from the initial con- 
ditions mentioned above that can be directly compared 
with the experimental data. Note that the procedure av- 
eraging over many samples without shift or rotation cor- 
responds to conventional initialization without the effect 
of eccentricity fluctuation. 



C. Monte Carlo Glauber model 

In the Monte Carlo version of the Glauber model, we 
first sample the positions of nucleons according to a nu- 
clear density distribution for two colliding nuclei. A 
nucleon-nucleon collision takes place if their distance d 
in the transverse plane (orthogonal to the beam axis) 
satisfies the condition 



d < 



'NN 



(14) 



where uj^^ ~ 41.94 mb is the inelastic nucleon-nucleon 
cross section at ^/snn ~ 200 GeV. The number of bi- 
nary collisions A^coii is obtained by counting total number 
of nucleon-nucleon collisions, and the number of partici- 
pants A^part is the total number of nucleons which collide 
at least once. 

In the ordinary MC-Glauber model, the following 
Woods-Saxon distribution is employed as a nuclear den- 
sity profile 



Pws 



exp I 



(15) 
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6.38 (4.20641) 



where po = 0.1695 (0.1686) fm 
fm d = 0.535 (0.5977) fm for a gold (copper) nucleus 
[iH . IH] . If one assumes a profile of nucleon as the delta 
function, the resulting nuclear density is nothing but the 
Woods-Saxon distribution 



p(f) = 6'^^^x - xo)pws{xo)d^xo. 



(16) 



However, in the case of finite size profile for nucleons as 
assumed in the MC-Glauber model, the nuclear density 
is no longer the same as the Woods-Saxon distribution. 

p{x) = J A{x - xo)pws{xo)d^xo, (17) 



A{x- xq) 
Vn 



(rN- I X - 



xo I) 



AT:r% 



Vn 

rN 



'NN 



(18) 
(19) 



This is illustrated in Fig. [TJ Obviously, nuclear surface 
is more diffused in both gold and copper nuclei due to 
the finite nucleon profile in Eq. (fT8|) . As a matter of fact, 
eccentricity becomes smaller by ~ 10% [lO] unless one 
adjusts the Woods-Saxon parameters according to the 
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D. Monte Carlo KLN model 

In the MC-KLN model, the saturation scale for a nu- 
cleus A at a transverse coordinate x±^ in each sample is 
given by 



,2 TAixj 



Ta,o 



(21) 



and similarly for a nucleus B. Here, parameters Qjg = 2 
GeV^, Ta,o = 1-53 frn'^, xq = 0.01, and A = 0.28 are 
used. Thickness function at each transverse coordinate 
is obtained by counting the number of wounded nucleons 
within a tube extending in the beam direction with 
radius r = \/<^nn/'^ from each grid point: 



FIG. 1: (Color online) Nuclear density as a function of nu- 
clear radius. Solid lines show nuclear density distribution for 
gold and copper nuclei in which a finite nucleon profile is im- 
plemented and positions of nucleons are sampled according 
to the Woods-Saxon distribution with default parameter sets. 
Dashed lines show the Woods-Saxon distribution with default 
parameter sets. 



finite nucleon profile. So we parameterize the distribu- 
tion of nucleon positions to reproduce the Woods-Saxon 
distribution with default parameters in Eq. (jlSp [2^ . In 
our MC-Glauber model, we find the default Woods-Saxon 
distribution is reproduced by a larger radius parameter 
ro = 6.42 (4.28) fm and a smaller diffuseness parame- 
ter d = 0.44 (0.50) fm (i.e., sharper boundary of a nu- 
cleus) for a gold (copper) nucleus than the default param- 
eters. This kind of effect exists almost all Monte Carlo 
approaches of the collisions including event generators. 
If one wants to discuss eccentricity and elliptic flow co- 
efficient V2 within ~ 10% accuracy, this effect should be 
taken into account. 

We assume that the initial entropy profile in the trans- 
verse plane at midrapidity is proportional to a linear com- 
bination of the number density of participants and that 
of binary collisions in the Glauber model 



dS 



Tadxdydrj, 
C_ fl 

To 



SdN, 



part 



.dN, 



coll 



d2 



d^ 



(20) 



Parameters C = 19.8 and 5 — 0.14 have been fixed 
through comparison with the centrality dependence of 
multiplicity data in Au-|-Au collisions at RHIC [1^ by 
pure hydrodynamic calculations with Tdec = 100 MeV. 
Notice that the total multiplicity of the produced parti- 
cle does not change during the late hadronic stage due 
to the chemical freezeout in the EOS. 



N 



'NN 



(22) 



For each generated configuration of nucleons in col- 
liding nuclei, the fcT-factorization formula is applied at 
each transverse coordinate to obtain the distribution of 
produced gluons locally. We apply the Kharzeev-Levin- 
Nardi (KLN) approach [i^l in the version previously em- 
ployed in [43|. In this approach, the distribution of gluons 
at each transverse coordinate Xj_ produced with rapidity 
y is given by the fcr-factorization formula [46t 



dNn 



2^ 

Cf 



d^x±dy Cf J Pt J 4 

X (t)Aixi,{pT+kTf/4;x^) 
X (t)Bix2,{pT-kTf/4;x±_), 



(23) 



where xi^2 =Pt exp(±y)/Y^ and pT is the transverse mo- 
mentum of the produced gluons. We choose an upper 
limit of lOGeV/c for the pT integration. For the uninte- 
grated gluon distribution function we use 



4>Aix, k'^;xj_) 



kCf (1 - .t)^ 



(24) 



where Cf = {N^ - l)/(2iVc). The parameter = 1.75 
is chosen for the overall normalization of the gluon mul- 
tiplicity in order to fit the multiplicity data in Au+Au 
coUisions at RHIC (isj . 

As an initial condition for hydrodynamical calcula- 
tions, the initial entropy density in the transverse plane 
is obtained by 



3.6 



dNa 



Tod^x±dr]s 



(25) 



y=Vs=0 



Here we identify gluon's momentum rapidity y with 
space-time rapidity rjs. 
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FIG. 2: (Color online) Centrality dependences of charged particle multiplicity in Au+Au (solid line) and Cu+Cu (dashed line) 
collisions at y^sjvjv ~ 200 GeV. (a) Initial conditions taken from the MC-Glauber model, (b) Initial conditions taken from the 
MC-KLN model. Experimental data are from Refs. [H, 




FIG. 3: (Color online) Eccentricity as a function of A'part in (a) Au+Au and (b) Cu+Cu collisions. Solid, dotted, dashed, and 
dash-dotted lines correspond to Epart from the MC-KLN model, Epart from the MC-Glauber model, erp from the MC-KLN 
model, and erp from the MC-Glauber model, respectively. 



III. RESULTS 

A. Centrality dependence of multiplicity 

Centrality dependences of multiplicity in the MC- 
Glauber model are compared with the PHOBOS data 
di, li^l in Fig. [2] (a). Centrahty dependence in Au+Au 
collisions is well reproduced with a two component (soft 
+ hard) model with a small fraction of hard compo- 
nent S = 0.14. The MC-KLN model with the setting 
mentioned in the previous section also gives a reason- 
able agreement with the data as shown in Fig. [2] (b). 
We did not tune initial parameters in Cu+Cu collisions. 
The multiplicity in the Glauber model depends only on 



.^part regardless of the collision system, while a non- 
trivial A'part dependence appears in the MC-KLN model: 
At a fixed A'part, {dN^h/dii)/{Np^,t/2) in Cu+Cu colli- 
sions is larger than the one in Au+Au collisions. All 
results are obtained by solving the hydrodynamic equa- 
tions below to Tdoc = 100 MeV. In the actual calculations 
of elliptic flow coefficients, we replace the hadronic flu- 
ids with the hadronic gases utilizing a hadronic cascade 
model. Nevertheless, the centrality dependence of the 
multiplicity is still within error bars (48| . 



6 



B. Eccentricity 

The eccentricities as functions of A^part with or with- 
out eccentricity fluctuations arc compared in the MC- 
Glauber and MC-KLN models in Au+Au (Fig. [3] (a)) 
and Cu+Cu (Fig. [3] (b)) collisions. The impact param- 
eter range, the number of participants, the eccentricity, 
and the transverse area are summarized in Appendix |^ 
for each centrality from 0-5% to 60-70% in Au-|-Au (Table 
m and Cu-hCu (Table |TJ) collisions at ^/IWn = 200 GeV. 
In semicentral Au-|-Au collisions (10-50% centrality), the 
effect of initial fluctuations enhances eccentricity param- 
eter by 8-11% (5-8%) in the MC-Glauber (MC-KLN) 
model. The enhancement factor epart/^RP is the largest 
at the very central bin (0-5% centrality): epart/sRP ~ 
1.83 in the MC-Glauber model and 1.53 in the MC-KLN 
model. A qualitatively similar behavior is observed in 
Cu-|-Cu collisions as shown in Fig. [3](b). However, it is 
quantitatively different: the enhancement factor is 1.51- 
1.74 (1.36-1.50) in the MC-Glauber (MC-KLN) model in 
semicentral collisions (10-50% centrality). In the very 
central events (0-5% centrality), the factor reaches 4.20 
(3.36) in the MC-Glauber (MC-KLN) model. Thus the 
resultant elliptic flow coefficient is expected to be en- 
hanced by almost the same amount of the factor due to 
eccentricity fluctuations especially in Cu+Cu collisions. 
We will demonstrate it by using a dynamical model in 
the next subsection. 

£rp in the MC-Glauber model is consistent with our 
previous result in which we used the standard eccentric- 
ity [l3| . Notice that the consistency is achieved only after 
retuning the Woods-Saxon parameters due to finite nu- 
cleon profiles discussed in the previous section, erp in 
the MC-KLN model is slightly smaller [2^ than the rc- 
suh from the naive KLN model In the MC-KLN 

model, there exists the minimum saturation scale which 
is nothing but the saturation scale of a nucleon. On the 
other hand, the saturation scale in the naive KLN model 
can be arbitrary small below to Aqcd, which leads to 
a sharper transverse profile and large eccentricity. We 
found eccentricity is reduced by ^ 15% in the MC-KLN 
model compared to the KLN model, which would affect 
the elliptic flow coefficients at the same amount. 



C. Elliptic flow coefficient 

Figure [4] shows the centrality dependence of V2 for 
charged hadrons at mid-rapidity (| 77 |< 1) in the Glauber 
model initialization in (a) Au-|-Au and (b) Cu-|-Cu col- 
lisions at y/SNN = 200 GeV. Experimental data [i^ are 
reasonably reproduced in Au-|-Au collisions. The effect 
of eccentricity fluctuations is not significant in Au+Au 
collisions, whereas V2 is largely enhanced in Cu+Cu col- 
lisions due to fluctuation effects. These tendencies are 
also expected from the results of initial eccentricity as 
shown in Fig. [S] Interestingly elliptic flow coefficients V2 
in the Glauber model initial conditions still slightly un- 



dershoot the experimental data, in particular, in Cu+Cu 
collisions even with fluctuation effects. 

FigureOis the same as Fig.|4]but the initial conditions 
are taken from the MC-KLN model. Again, the effect 
of eccentricity fluctuations is small in Au+Au collisions 
but is large in Cu+Cu collisions. Due to larger initial ec- 
centricity in the MC-KLN model than the MC-Glauber 
model, the results are somewhat larger than the exper- 
imental data in peripheral Au+Au collisions. Whereas, 
we reasonably reproduce the V2 data in Cu+Cu collisions. 
It should be noted that the results in Au+Au collisions 
lie systematically below the ones in Fig. 2 of Ref. [l3| 
which employed the naive KLN model. This is due to the 
reduction of eccentricity by employing the fKLN model. 

In both cases, the results using a more sophisticated 
EOS such as the one from the numerical simulations of 
the lattice QCD ^ would be demanded. The EOS 
from the lattice simulations has a cross-over behavior 
rather than the flrst order phase transition. Therefore 
it is harder in the vicinity of pseudo-critical temperature 
Tc ~ 190 MeV and softer in the high temperature region 
up to T ^ 2Tc than the one employed in the present 
study. Because of the above reasons, whether flnal V2 
is enhanced or reduced in comparison with the current 
results would be non-trivial in the case of the EOS from 
the lattice QCD. Detailed studies on the dependences of 
V2 on the EOS will be reported elsewhere. 



IV. CONCLUSIONS 

We calculated the elliptic flow parameter as a function 
of the number of participants in the QGP hydro plus 
the hadronic cascade model and found that the effect of 
eccentricity fluctuations is visible in very central and pe- 
ripheral Au+Au collisions and is quite large in Cu+Cu 
collisions. This strongly suggests that the effect of ec- 
centricity fluctuations is an important factor which has 
to be included in the dynamical model for understand- 
ing of the elliptic flow data and for precise extraction of 
transport properties from the data. 

We also found that a flnite nucleon size assumed in the 
conventional Monte Carlo approaches reduces initial ec- 
centricity by ^10% with a default Woods-Saxon param- 
eter set. This requires reparametrization of the nuclear 
radius and the diffuseness parameter to obtain the ac- 
tual nuclear distribution in the case of the finite profile 
of nuclcons as implemented in the Monte Carlo Glauber 
model. 

In the case of the Glauber type initialization, the re- 
sults still undershot the experimental data a little in both 
Au+Au and Cu+Cu collisions even after inclusion of ec- 
centricity fluctuation effects. If one implemented the vis- 
cosity in the QGP phase, the results of V2 would get 
reduced. So there is almost no room for the viscosity in 
the QGP stage to play a role in the Glauber initial condi- 
tions within the hybrid approach with the ideal gas EOS 
in the QGP phase. 
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FIG. 4: (Color online) Centrality dependences of V2 for charged particles at mid-rapidity in (a) Au+Au and (b) Cu+Cu 
collisions at y^ijvjv ~ 200 GeV with the MC-Glauber model initial conditions are compared with PHOBOS data (filled plots) 
[49I ]. Open circles (squares) are results with (without) eccentricity fluctuation. 
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FIG. 5: (Color online) Centrality dependences of V2 for charged particles at mid-rapidity in (a) Au-f Au and (b) Cu-|-Cu 
collisions at ^snn = 200 GeV with the MC-KLN model initial conditions are compared with PHOBOS data (filled plots) [igj ]. 
Open circles (squares) are results with (without) eccentricity fluctuation. 



On the other hand, we overpredicted V2 in peripheral 
(^'^part < 150) Au4-Au colhsions in the CGC model. Vis- 
cous effects in the QGP phase could reduce the V2 and en- 
able us to reproduce the data in Au-|-Au collisions in this 
case. However, the results are already comparable with 
the data in Cu-|-Cu collisions even though the number 
of participants is almost the same as that in peripheral 
Au-|-Au collisions. So it would be non-trivial whether the 
same viscous effects also give the right amount of V2 in 
Cu4-Cu collisions. 

So far, one has been focusing on comparison of hydro- 
dynamic results with V2 data only in Au-f-Au collisions. 
The experimental data in Cu-f Cu collisions also have a 
strong power to constrain the dynamical models. There- 



fore, simultaneous analysis of V2 data in both Au+Au 
and Cu-f Cu collisions will be called for in future hydro- 
dynamic studies. 
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For a given range of centrahty Cmin < c < Cmax? the max- 
imum and minimum impact parameters can be defined 
from Eq. jAl]) as c(6niax) = c^ax and c(6i„in) = Cmin, 
respectively. In generating initial profiles, an impact pa- 
rameter is randomly chosen in b^i^ < b < 6niax for each 
centrality bin with a probability 



APPENDIX A: CENTRALITY DEPENDENCE OF 
ECCENTRICITY AND TRANSVERSE AREA 



P{b)db = 



2b 



62 - 62 . 

"max '^mm 



db. 



(A3) 



We obtain centrality in our model as follows. Using 
the Glauber model, one calculates centrality, namely the 
fraction of inelastic cross section as a function of impact 
parameter, 

cib) = h — — (Al) 



2^/„°°6'd6'(l-e— ««^^^(^')) 



TAA{b) = J d^x±TA (^x - ^,y^ Ta 



(A2) 



Impact parameter ranges, the resultant number of par- 
ticipants, eccentricity with respect to the reaction plane, 
participant eccentricity, and transverse area are summa- 
rized in Table [T] for Au+Au collisions and Table |TT] for 
Cu+Cu collisions. These parameters would be very use- 
ful to sec whether data reach the so-called hydrody- 
namic limit or not |5l| . It might be possible to divide 
all events into centralitics according to the Monte Carlo 
results of, e.g., the iVpart distribution. 
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TABLE I: Centrality dependence of eccentricity and transverse area in Au+Au collisions. 



Centrality(%) 


0-5 


5-10 


10-15 


15-20 


20-30 


30-40 


40-50 


50-60 


60-70 


6min (fm) 


0.0 


3.3 


4.7 


5.8 


6.7 


8.2 


9.4 


10.6 


11.6 


&max (fm) 


3.3 


4.7 


5.8 


6.7 


8.2 


9.4 


10.6 


11.6 


12.5 


A^part 


352 


295 


245 


204 


154 


104 


65.1 


36.8 


18.8 


MC-Glauber 
^RP 


0.0446 


0.120 


0.183 


0.233 


0.292 


0.348 


0.389 


0.405 


0.398 


MC-Glauber 
part 


0.0818 


0.145 


0.204 


0.254 


0.316 


0.380 


0.433 


0.473 


0.497 


^.MC-Glaubcr (f^2) 


23.4 


20.5 


18.0 


16.0 


13.5 


10.9 


8.69 


6.78 


5.07 


^MC-Gla.bcr (f^2) 


23.4 


20.5 


18.0 


16.0 


13.5 


10.9 


8.65 


6.73 


5.05 


MC-KLN 
^RP 


0.0671 


0.165 


0.239 


0.298 


0.367 


0.435 


0.485 


0.509 


0.506 


MC-KLN 
part 


0.103 


0.185 


0.257 


0.316 


0.387 


0.461 


0.522 


0.571 


0.599 


gMC-KLN (f^2) 


23.8 


20.1 


17.3 


15.1 


12.4 


9.67 


7.42 


5.55 


3.96 


5-^-"^" (W) 


23.7 


20.1 


17.2 


15.0 


12.3 


9.55 


7.26 


5.33 


3.67 



TABLE II: Centrality dependence of eccentricity and transverse area in Cu+Cu collisions. 



Centrality(%) 
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MC-Glaubcr 
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0.257 
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MC-Glaubor 
part 


0.115 
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0.434 
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oMC-Glaubcr fr 2\ 

^RP llm ) 


12.1 
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10.1 


9.15 


7.98 
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3.00 
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1.89 



10 



[27] C. Nonaka and S. A. Bass, Nucl. Phys. A 774, 873 (2006); 

Phys. Rev. C 75, 014902 (2007). 
[28] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, and 

H. Stocker, Phys Rev. C 78, 0449 01 (2008); H. Petersen 

and M. Bleicher. larXiv:0901.382T1 [nucl-th]. 
[29] K. Werner, T. Hirano, lu Karpenko, T. Pierog, S. Por- 

teboeuf, M. Bleicher, and S. Haussler, J. Phys. G: 

Nucl. Part. Phys. (in press). 
[30] F. Cooper and G. Frye, Phys. Rev. D 10, 186 (1974). 
[31] I. G. Bearden et al, [BRAHMS Collaboration], Phys. 

Rev. Lett. 93, 102301 (2004). 
[32] T. Hirano, Phys. Rev. C 65, 011901 (2002). 
[33] T. Hirano and K. Tsuda, Phys. Rev. C 66, 054905 (2002). 
[34] H. Bebie, P. Gerber, J. L. Goity, and H. Leutwyler, Nucl. 

Phys. B 378, 95 (1992). 
[35] N. Arbex, F. Grassi, Y. Hama, and O. Socolowski Jr., 

Phys. Rev. C 64, 064906 (2001); W. L. Qian, R. Andrade, 

F. Grassi, O. Socolowski Jr., T. Kodama, and Y. Hama, 

Int. J. Mod. Phys. E 16, 1877 (2007). 
[36] D. Teaney, nucl-th/0204023 

[37] P. F. Kolb and R. Rapp, Phys. Rev. C 67, 044903 (2003). 

[38] P. Huovinen, a rXrvT0710. 4379 [nucl-th]. 

[39] Y. Nara, N. Otuka, A. Ohnishi, K. Niita, and S. Chiba, 

Phys. Rev. C 61, 024901 (2000). 
[40] R. P. G. Andrade, F. Grassi, Y. Hama, T. Kodama, and 

W. L. Qian, Phys. Rev. Lett. 101, 112301 (2008). 



[41] M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, 
Ann. Rev. Nucl. Part. Sci. 57, 205 (2007). 

[42] D. Kharzeev and M. Nardi, Phys. Lett. B 507, 121 
(2001); D. Kharzeev and E. Levin, Phys. Lett. B 523, 
79 (2001); D. Kharzeev, E. Levin, and M. Nardi, Phys. 
Rev. C 71, 054903 (2005); D. Kharzeev, E. Levin, and 
M. Nardi, Nucl. Phys. A 730, 448 (2004). 

[43] T. Hirano and Y. Nara, Nucl. Phys. A 743, 305 (2004). 

[44] H. De Vries, C. W. De Jager, and C. De Vries, Atom. 
Data Nucl. Data Tabl. 36, 495 (1987). 

[45] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. 
C 65, 061901 (2002). 

[46] L. V. Gribov, E. M. Levin, and M. G. Ryskin, Phys. 
Rept. 100, 1 (1983). 

[47] R Alver et al. [PHOBOS Collaboration], 

larXiv:0808.1895 [nucl-ex]. 

[48] T. Hirano, \]. Heinz, D. Kharzeev, R. Lacey, and Y. Nara, 
J. Phys. G: Nucl. Part. Phys. 34, S879 (2007). 

[49] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. 
C 72, 051901 (2005); B. Alver et al. [PHOBOS Collabo- 
ration], Phys. Rev. Lett. 98, 242302 (2007). 

[50] M. Cheng et al, Phys. Rev. D 77, 014511 (2008); 
A. Bazavov et a/.. .arXiv:0 903.4379 [hep-lat]. 

[51] H. J. Drescher, A. Dumitru, C. Gombeaud, and J. Y. Ol- 
litrauh, Phys. Rev. C 76, 024905 (2007). 



